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Numerical  modelling  of  redox  flow  battery  (RFB)  systems  for  energy  storage  applications  allows  the 
technical  performance  of  different  designs  to  be  predicted  without  costly  lab,  pilot  and  full-scale  testing. 
A  one-dimensional  numerical  model  has  been  developed  for  RFB  systems  with  bipolar  flow-by  electrodes, 
soluble  redox  couples,  and  recirculating  batch  operation.  Overpotential  losses  were  estimated  from  the 
Butler-Volmer  equation,  accounting  for  mass  transfer.  The  effects  of  cross-membrane  solvent  transport 
and  self-discharge  were  also  considered.  The  model  predicted  the  variation  in  concentration  and  current 
along  the  electrode  and  determined  the  charge-discharge  efficiency,  energy  density  and  power  density. 
The  model  was  validated  using  data  obtained  from  a  pilot-scale  bromide-polysulphide  (PSB)  system 
commercialised  by  Regenesys  Technologies  (UK)  Ltd.  Electrochemical  rate  constants  were  obtained  by 
fitting  the  model  results  to  the  experimental  data,  and  values  of  4  x  10-7  and  3  x  10-8  ms-1  were  found 
for  the  bromide  and  sulphide  electrolytes  on  the  activated  carbon  electrodes.  The  model  was  able  to 
predict  cell  performance,  species  concentration,  current  distribution  and  electrolyte  deterioration  for  the 
Regenesys  system. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  demand  for  utility-scale  energy  storage  is  growing  due  to  an 
increasing  need  for  load-smoothing  to  allow  the  decommissioning 
of  old  power  plant,  to  improve  power  quality  for  the  growing 
digital  economy,  to  improve  security  of  supply  through  black-start 
capability  and  to  smooth  the  output  from  renewable  energy 
technologies.  Redox  flow  batteries  (RFBs)  have  been  investigated 
for  many  years  as  chemical  stores  of  electrical  energy  [1,2],  and  are 
the  closest  storage  technology  to  widespread  commercialisation 
[3].  RFBs  have  numerous  advantages  over  other  batteries,  including 
a  separation  of  the  energy  and  power  rating,  modular  systems, 
repeatable  cyclic  behaviour  and  the  use  of  benign  chemicals  at 
atmospheric  temperature  and  pressure.  Redox  couples  currently 
under  development  for  use  in  RFBs  include  polysulphide-bromine 
(PSB)  [4-7],  vanadium-vanadium  [8,9],  vanadium-polyhalide  [10], 
cerium-zinc  [11-13]  and  lead-lead  [2].  Numerical  modelling  of 
redox  flow  battery  systems  for  energy  storage  applications  allows 
the  technical  and  commercial  performance  of  different  designs  to 
be  predicted  without  costly  lab,  pilot  and  full-scale  testing.  Limited 
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modelling  of  RFB  systems  has  been  carried  out  to  date.  A  study  of 
the  iron-chromium  system  included  reaction  kinetics,  mass  trans¬ 
fer,  ohmic  losses  and  hydrogen  evolution,  and  found  an  optimum 
electrode  thickness  and  electrolyte  flow  rate  [14].  A  model  of  the 
zinc-ferricyanide  couple  considered  only  mass-transport-limited 
conditions  [15].  An  investigation  of  the  bromide  half-couple  found 
that  a  couple  of  reaction  mechanisms  could  account  for  observed 
behaviour  [16].  A  review  of  modelling  work  on  the  zinc-bromide 
system  has  reported  recommendations  for  improving  cell  design 
[17]. 

In  this  paper  we  develop  a  numerical  model  of  a  RFB  system,  and 
apply  this  to  the  Regenesys  Technologies  Ltd.  pilot-scale  PSB-based 
RFB.  The  PSB  RFB  has  not  previously  been  studied  by  numerical 
simulation.  Numerical  modelling  can  be  used  to  obtain  key  param¬ 
eters  such  as  the  electrochemical  rate  constants  for  the  reactions. 
Furthermore,  once  validated,  the  model  can  be  used  to  evaluate  and 
optimise  the  design  and  performance  of  a  full-scale  commercial  RFB 
system.  This  is  the  subject  of  a  following  paper  [18]. 

2.  Method 

2.1.  Cell  layout 

The  model  centred  on  a  single,  vertical,  membrane-divided 
regenerative  electrochemical  cell  as  shown  in  Fig.  1 .  The  electrolytes 
entered  at  the  bottom  of  the  cell  (to  allow  easy  filling,  promote 
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Nomenclature 

Nomenclature 

A  electrode  area  (m2  cell-1 ) 

/Ichan  channel  flow  area  (m2  cell-1 ) 
de  electrode  hydraulic  diameter  (m) 

D  diffusion  coefficient  (m2  s-1 ) 

E  voltage  (V) 

E°  standard  potential  (V) 

F  Faraday’s  constant  (96,485  C  mol-1 ) 

H  number  of  segments 

I  current  (A) 

i  current  density  (A  m-2 ) 

i0  exchange  current  density  (A  m-2 ) 

id  dimensionless  current  i/k 

iL  limiting  current  density  (A  m-2 ) 

iir  limiting  current  ratio 

ir  current  ratio 

J  species  flux  across  the  membrane  (mol  m-2  s-1 ) 

km  mass  transfer  coefficient  (m  s-1 ) 

ks  kinetic  rate  constant  (m  s-1 ) 

N  number  of  cells 

Njt  number  of  iterations 

0  oxidised  species* 

[P]  concentration  of  species  P  (mol  dm-3 ) 

R  reduced  species* 

R  universal  gas  constant  (8.314J  mol-1 1<-1 ) 

R  resistance  (£2) 

l  thickness  (m) 

t  time  (s) 

A  t  time  increment  ( s ) 

T  temperature  (I<) 

Tcomp  computational  time  (s) 

v  electrolyte  velocity  (m  s-1 ) 

V  electrolyte  volume  (m3) 

z  number  of  electrons  transferred  per  occurrence  of 

overall  reaction* 

Zp  stoichiometric  coefficients  of  species  P 

Greek  symbols 

a  charge  transfer  coefficient 

s  relative  error 

r]  overpotential  (V) 

r]f  dimensionless  overpotential 

k  conductivity  ( S  m- 1 ) 

Ah  convergence  criterion  of  segment  iteration  (A  m-2 ) 

Ae  convergence  criterion  of  electrode  iteration  (V) 

v  kinematic  viscosity  (m2  s-1 ) 

a  state  of  charge 

Subscripts 

a  anodic  value 

A  half-couple  A 

ave  average  value 

B  half-couple  B 

Br  value  for  the  bromide  electrolyte 

bulk  value  in  electrolyte  bulk 

c  cathodic  value 

CELL  cell  value 

ch  value  on  charge  cycle 

dis  value  on  discharge  cycle 

entry  value  at  electrode  entry 

eqm  equilibrium  value 

exit  value  at  electrode  exit 


Est  estimated  value 

f  value  for  forward  reaction 

h  segment  value 

M  membrane 

P  value  for  general  species  P 

S  value  for  the  sulphide  electrolyte 

T  tank  value 

Dimensionless  numbers 
Re  Reynolds  number  ( vde/v ) 

Sc  Schmidt  number  (v/D) 

Sh  Sherwood  number  (kmde/D) 


good  mixing  and  even  conditions  across  the  cell,  and  allow  gases  to 
escape  easily),  were  pumped  upwards  where  the  dissolved  active 
species  reacted  at  flow-by  electrodes,  and  left  at  the  top  to  be  stored 
in  tanks.  On  charge,  redox  couple  A  was  oxidised  at  the  anode  and 
couple  B  was  reduced  at  the  cathode  while  cations  crossed  the 
membrane  to  maintain  electro-neutrality,  with  the  reverse  occur¬ 
ring  on  discharge.  Multiple  cells  were  connected  together  in  parallel 
using  bipolar  electrodes  to  form  stacks,  and  electrolytes  were  fed 
into  and  collected  from  the  cell  using  manifolds.  The  same  tank  was 
used  to  store  the  reactants  and  products,  so  that  the  electrolytes 
recirculated  with  the  concentration  of  product  steadily  increasing 
over  time.  A  constant  current  was  supplied,  and  the  cell  voltage 
changed  during  the  charging  cycle  as  the  concentrations  of  reac¬ 
tants  and  products  varied. 

2.2.  Model  scope  and  assumptions 

The  RFB  numerical  model  developed  in  this  paper  performed 
the  following  functions: 

•  Evaluation  of  mass  transport  (from  the  electrolyte  bulk  to  the 
electrode  surface)  and  reaction  kinetics  (described  by  the  Butler- 
Volmer  equation  with  lst-order  reaction  rates  and  equal  rate 
constants  for  anodic  and  cathodic  processes)  as  rate-determining 
processes. 

•  Estimation  of  the  variation  in  concentration,  overpotential,  cur¬ 
rent  density,  exchange  current  density  and  limiting  current 
density  up  the  electrode  in  a  one-dimensional  model. 

•  Calculation  of  the  variation  in  cell  performance  during 
charge-discharge  cycles  and  overall  system  characteristics 
including  energy  efficiency,  power  density  and  energy  density. 

•  Consideration  of  different  operating  conditions  and  electrolyte 
systems,  e.g.  variable  redox  couple,  applied  current  density,  power 
rating,  operating  temperature,  catalyst,  cycle  length,  species  con¬ 
centration,  electrolyte  velocity,  electrode  area,  stack  size,  tank 
volume,  electrolyte  conductivity  and  membrane  conductivity. 

•  Estimation  of  the  effect  of  the  diffusion  of  electro-active  species 
and  electro-osmotic  flow  of  water  across  the  membrane. 

The  main  assumptions  used  in  the  model  are  listed  below: 

•  Single-step  electrochemical  reactions  involving  dissolved 
electro-active  species  were  assumed  to  occur  at  the  electrodes. 

•  Electrochemical  kinetics  were  assumed  to  be  described  by  the 
Butler-Volmer  equation,  with  the  anodic  and  cathodic  charge 
transfer  coefficients  summing  to  one. 

•  Electrochemical  rate  constants  and  mass  transfer  coefficients 
were  assumed  to  be  approximately  constant,  i.e.  the  effect  of 
temperature,  electrode-fouling,  etc.  was  not  included. 
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Fig.  1.  Cell  layout  and  electrochemical  conventions. 


•  The  effects  of  electro-migration  and  mass  transport  away  from 
the  electrode  were  assumed  to  be  negligible. 

•  The  effects  of  adsorption,  capacitance,  electrodes/endplates  resis¬ 
tivity  and  shunt  currents  were  assumed  to  be  negligible. 

•  The  current  efficiency  was  assumed  to  be  100%,  i.e.  side  reactions 
were  not  considered. 

•  Plug  flow  conditions  were  assumed  to  occur  in  the  cells,  with  a 
constant  electrolyte  velocity  in  both  half-cells  and  during  cycles. 

•  Conditions  were  assumed  to  be  the  same  in  each  cell  in  a  stack, 
so  that  the  performance  of  a  stack  of  cells  could  be  determined 
directly  from  the  performance  of  a  single  cell. 

2.3.  Model  equations 

The  model  assumed  that  simple,  single-step,  single-species  elec¬ 
trochemical  reactions  occurred  at  each  electrode: 

z00  +  ze~  ^zrR  (1) 

The  cell  voltage  was  calculated  from  equilibrium  potentials  at 

the  electrodes,  accounting  for  the  main  losses  in  the  cell.  During 

charging  the  cell  voltage  was  calculated  as: 

^CELL  =  £eqm,A  “  ^eqm.B  +  dA  ~  ?7b  +  (f^A  +  ^B  +  IRm)  (2) 

On  discharging,  the  cell  voltage  was  calculated  as: 


resulting  at  each  electrode  was  determined  by  inverting  the  Butler- 
Volmer  equation,  accounting  for  the  effect  of  mass  transport  on  the 
active  species  concentration  at  the  electrode: 

i= _ exp{(aazF/RT)q}~ exp  {-{aczF/RT)q} _ 

(1/t'o)  +  ( 1  / iLa ) exp{(o'azf/ RT)rj}  -  {J/ilc)exp{-{aczF/RT)q} 

where  the  exchange  current  density  z0,  anodic  limiting  current  den¬ 
sity  iLa,  and  cathodic  limiting  current  density  zLc  were  defined  as: 

io=zFks[0\ab^k[R\ah^k  (7) 

*La  —  ^^m,R[^]bulk  (8) 

*Lc  —  — zF/cm,o[G]bUik  (9) 

In  Eqs.  (6)-(9),  both  z0  and  zR  equal  1.  No  adequate  rate  expres¬ 
sion  is  currently  known  for  electrochemical  reactions  with  values 
for  z0  and  zR  other  than  1,  and  it  has  been  observed  that  com¬ 
plex  reaction  mechanisms  often  show  simple  behaviour  and  can 
be  treated  by  relations  of  this  sort  [19].  It  is  likely  that  values  of  Zo 
and  zR  other  than  1  will  affect  the  scaling  but  not  the  form  of  the 
Butler-Volmer  equation,  and  hence  it  is  reasonable  to  neglect  their 
influence  on  Eqs.  (6)-(9). 

The  inversion  of  Eq.  (6)  was  achieved  using  approximations  to 
the  Butler-Volmer  equation  rather  than  iteration  in  order  to  reduce 
computational  time.  Eq.  (6)  was  written  in  dimensionless  format: 


^CELL  =  £eqm,A  “  ^eqm.B  +  dA  ~  dB  ~  +  f^B  +  IRm)  (3) 


The  equilibrium  potentials  at  each  electrode,  Eeqm,A  and  Eeqm,B 
were  derived  from  the  standard  potential  and  species  concentra¬ 
tions  (assuming  unit  activity  coefficients)  as  given  in  the  Nernst 
equation: 


E  —  E°  RT  In  [RfR 

teqm  -  b  ~  W  [Op 


(4) 


The  electrolyte  and  membrane  resistances  RA,  RB  and  RM  were 
derived  from: 


(5) 


where  l  is  the  thickness  (membrane  or  electrolyte),  /c  is  the  conduc¬ 
tivity  and  A  is  the  active  area  of  the  electrode.  The  overpotentials  q 


=  exp{o:f??f)  -  exp{-(l  -  Of  )qf) 

*Lf  «rf  +  exp{o!f?7f}  +  iLr  exp{-(l  -  af)r/fj 


(10) 


The  forward  process  was  either  anodic  or  cathodic  depending 
on  the  direction  of  the  net  reaction.  The  dimensionless  current  idf 
varied  between  0  and  +1  and  gave  a  simple  indication  of  the  near¬ 
ness  of  the  current  to  mass-transport  limits.  Other  dimensionless 
quantities  were  defined  as: 


(11) 


(12) 

(13) 
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The  current  ratio  (ir f)  gave  an  indication  of  the  relative  rates 
of  mass  transport  and  reaction  kinetics.  Eqs.  (10)  and  (11)  define 
the  dimensionless  overpotential  as  a  function  of  four  positive 
quantities:  zdf,  af,  zrf  and  iiT. 

Three  approximations  were  used  to  invert  Eq.  (10):  the  Linear, 
mass-transport-limited  (MTL)  Tafel  and  Sigmoidal  approximations: 


df  =  irfidf 
1  , 

m=  —  In 

a  f 


jrfjdf 
.1  -  *df- 


*lf  =  ln 


+_ZLrjdf~ 

.  1  -  *df  . 


(14) 

(15) 

(16) 


The  Linear  approximation  was  used  for  conditions  of  low  zdf,  the 
MTL  Tafel  approximation  in  conditions  of  high  zdf,  and  the  Sigmoidal 
approximation  in  conditions  of  low  irf.  Further  details  of  the  inver¬ 
sion  method  are  described  elsewhere  [20],  where  it  has  been  shown 
that  the  error  in  overpotential  obtained  from  using  these  approxi¬ 
mations  was  of  the  order  of  10-5.  At  the  transition  points  between 
the  approximations  the  error  could  rise  to  around  10%,  although  the 
resulting  error  in  derived  quantities  such  as  cell  voltage  and  time- 
averaged  quantities  such  as  current  efficiency  would  be  much  less 
than  this. 


2.4.  Model  methodology 

The  model  methodology  is  shown  in  Fig.  2.  First,  the  cell  volt¬ 
age  FEst  (assumed  to  be  constant  up  the  electrode  length)  was 
estimated.  A  one-dimensional  numerical  model  of  the  current 
distribution  along  the  electrode  was  calculated  by  dividing  the 
electrode  into  H  segments  and  calculating  the  current  in  each 
segment  iteratively.  The  current  zh  in  an  individual  segment  was 
estimated,  and  used  to  find  the  species  exit  concentrations.  In  each 
segment  the  species  concentrations,  equilibrium  potentials  Feqm,A 
and  Eeqm.B,  overpotentials  r]A  and  z/B  and  the  ohmic  losses  in  the 
electrolyte  and  membrane  were  calculated.  The  segment  voltage  Fh 
was  calculated  from  Eqs.  (2)  and  (3),  and  an  interpolation-bisection 
iteration  of  zh  was  carried  out  until  the  discrepancy  between  Fh 
and  FEst  fell  below  a  segment  convergence  criterion  Ah.  This  itera¬ 
tion  was  repeated  for  all  the  segments  along  the  electrode,  and  the 
average  current  density  iave  found  from  the  resulting  current  dis¬ 
tribution.  A  second  iteration  was  performed,  varying  FEst  until  the 
discrepancy  between  iave  and  zCEll  was  less  than  an  electrode  con¬ 
vergence  criterion  Xe.  This  gave  the  cell  voltage  £CEll  at  time  t.  The 
above  calculation  process  was  repeated  for  the  length  of  the  charg¬ 
ing  cycle  and  for  subsequent  charging  and  discharging  cycles  using 
Faraday’s  law  and  a  mass  balance  to  calculate  the  concentration  of 
the  electrolytes  in  the  storage  tanks. 

Both  of  the  iteration  loops  in  the  model  were  executed  using 
an  interpolation-bisection  method.  Bisection  methods  converge  as 
long  as  initial  limits  are  chosen  either  side  of  the  solution,  and  as  the 
voltage-current  relationship  in  the  mixed  kinetic-limited/mass- 
transport-limited  region  of  the  Butler-Volmer  equation  frequently 
encountered  in  this  study  is  approximately  linear,  linear  interpola¬ 
tion  was  found  to  converge  quickly. 

3.  The  bromide-polysulphide  couple 

For  long-term  storage,  high  solubility,  low  toxicity  and  low  cost 
salts  are  needed  and  alkali  metal  salts  of  bromine,  chlorine  and  sul¬ 
phur  are  ideal.  These  salts  are  attractive  due  to  their  abundance 
and  availability  at  the  necessary  degree  of  purity  at  moderate  cost. 
The  sodium  bromide-polysulphide  redox  couple  is  attractive  for 
flow  batteries  because  it  presents  no  adverse  hazards  in  handling 
or  storage  in  its  uncharged  state,  whereas  other  redox  couples  are 
not  as  suitable  for  reasons  of  safety,  efficiency  and  cost  [5].  Bromine 


Fig.  2.  Flowchart  of  numerical  redox  flow  battery  model. 


has  long  been  investigated  as  an  electrolyte  for  redox  flow  batter¬ 
ies,  normally  with  zinc  as  the  other  electro-active  species  [21,22], 
and  the  PSB  couple  avoids  the  problem  of  handling  deposited  solid 
metal.  The  high  stability  of  the  bromide  ion  gives  it  a  high  poten¬ 
tial,  making  it  attractive  for  redox  flow  batteries.  Bromine  is  only 
sparingly  soluble  in  water  but  does  complex  as  soluble  tribromide 
ions  Br3_,  and  this  increases  its  solubility  at  room  temperatures. 
Polysulphide  has  not  been  considered  as  an  electrolyte  as  exten¬ 
sively  as  bromide.  The  dissolved  sulphur  appears  as  polysulphide 
ions  with  a  maximum  chain  length  of  five  sulphur  atoms,  and  usu¬ 
ally  as  S42-  or  S52-  ions.  The  maximum  solubility  of  I<2S  is  8.8  M 
I<2S  at  25  °C,  which  corresponds  to  2.2  M  I<2S4  when  discharged 
[23].  This  ability  to  store  large  amounts  of  dissolved  zero-valent 
sulphur  as  polysulphides  means  that  sulphur  solutions  are  good 
catholytes  [24].  The  standard  potential  for  the  bromide  and  poly¬ 
sulphide  reactions  are  +1.087  and  -0.508  V  respectively,  giving  an 
overall  standard  potential  of  1.59  V.  This  is  close  to  the  maximum 
achievable  in  water-based  systems  of  about  2  V,  but  means  that  gas 
evolution  occurs  at  both  high  and  low  states  of  charge.  The  state  of 
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charge  of  each  electrolyte  was  defined  as: 


4.1.  Stack  design 


[Br2] 

[Br2]  +  i[Br-] 

(17) 

[s2-]  -  1[S] 

[S2-]  +  [S] 

(18) 

Here  [S]  refers  to  the  total  concentration  of  zero- valent  sulphur. 
These  expressions  are  independent  of  electrolyte  volume.  The  fresh 
electrolytes  consisted  of  3.8  m3  of  4.5  M  NaBr  and  3.2  m3  of  1.0  M 
Na2S4.8,  with  VBr  greater  than  Vs  to  allow  for  cross-membrane  water 
transport.  It  can  be  simply  shown  that,  for  these  fresh  compositions, 
<rBr  and  crs  are  linked  by  [20]: 


75  1 

128aBr  +  96 


(19) 


The  conductivities  and  viscosities  of  the  electrolytes  have 
been  determined  at  35  °C  and  at  mid-charge:  /cBr  =  23  Sm-1, 
ks  =  29  S  nrr1 ,  vBr  =  0.76  x  10-6  m2  s-1  and  vs  =  1.18  x  10-6  m2  s_1 . 

A  number  of  mechanisms  have  been  proposed  for  the 
bromide-polysulphide  system  under  a  range  of  conditions  [25-31  ] 
including  tribromide  and  polysulphide  complexes,  supersulphide 
radicals  and  intermediary  steps.  In  this  study  a  simple  reaction 
mechanism  is  assumed: 


Br2(aq)  +2e  ^  2Br 

(20) 

S  +  2e"  ^  S2- 

(21) 

The  kinetic  data  required  for  the  model  included  rate  constants 
and  charge  transfer  coefficients.  For  the  bromide  couple,  the  rate 
constant  was  estimated  from  values  of  i0  [20]  to  be  10_4-10-3  ms-1 
on  platinum  electrodes  [25,27]  and  10_8-10~4 ms-1  on  vitreous 
carbon  electrodes  [26,32].  No  data  is  currently  available  for  the 
bromide  couple  alone  on  activated  carbon  electrodes,  but  this  is 
likely  to  be  similar  to  vitreous  carbon,  i.e.  a  value  in  the  range 
10_6-10-8  ms-1  is  reasonable.  These  authors  also  found  a  range 
of  charge  transfer  coefficients  from  Tafel  slopes,  and  it  seems  rea¬ 
sonable  to  use  values  for  both  aa  and  ac  of  0.5. 

For  the  sulphide/polysulphide  couple,  the  rate  constant  on  a 
range  of  metals  including  gold,  tungsten,  nickel,  graphite,  platinum 
and  cobalt  [29,33]  has  been  estimated  from  i0  values  to  vary  in  the 
range  10_1°-10-7  m  s-1  [20].  On  activated  carbon  [34],  the  rate  con¬ 
stant  has  been  estimated  to  be  10-7  m  s-1  [20].  Other  studies  found 
that  the  cathodic  reaction  had  a  high  overpotential  on  gold  [35],  that 
Qfa  on  graphite  was  0.5  [36],  and  that  the  rate  constant  is  lst-order 
[37].  A  range  of  charge  transfer  coefficients  have  been  estimated, 
and  it  seems  reasonable  to  use  values  for  both  aa  and  ac  of  0.5. 

A  study  of  the  PSB  system  as  a  whole  on  activated  carbon  found 
that  the  total  overpotential  for  a  PSB  system  at  40mAnrr2,  299  K, 
50%  charged  4.0  M  NaBr  and  1.3  M  Na2S4  electrolytes  is  191  mV  [7]. 
If  half  of  this  is  assumed  to  be  due  to  the  membrane  and  electrolyte 
resistive  losses  [38],  if  the  overpotential  is  assumed  to  caused  by 
activation  only,  and  if  the  bromide  rate  constant  is  assumed  to 
be  an  order  of  magnitude  greater  than  the  sulphide  rate  constant 
(see  Section  4),  then  ks  Br  and  /<s  S  can  be  estimated  as  5  x  10-7  and 
5  x  10-8  ms-1  respectively. 


4.  Modelling  the  pilot-scale  bromide-polysulphide  system 

A  pilot-scale  PSB  RFB  system  was  tested  by  Regenesys  Technolo¬ 
gies  Ltd.,  and  data  from  this  system  was  used  to  validate  the  model. 
This  system  used  a  stack  of  200  bipolar  electrochemical  cells,  called 
the  XL200. 


The  XL  cell  considered  in  this  paper  had  an  electrode  width 
and  length  of  66.8  and  108.0  cm  respectively  with  a  total  pro¬ 
jected  area  of  0.67  m2.  The  electrode-membrane  gap  dg  averaged 
about  0.95  mm  to  give  a  flow  cross-sectional  area  of  6.5  cm2. 
The  bipolar  electrodes  consisted  of  an  electrically  conductive, 
non-porous  high-density  polyethylene  (HDPE)/graphite  compos¬ 
ite  core  on  which  particles  of  mixed  polyvinylidene  difluoride 
(PVDF)  and  activated  carbon  were  compressed  on  both  sides  to 
increase  the  active  surface  area  [39].  These  electrodes  were  laser- 
welded  into  an  HDPE  frame  and  were  about  5.5  mm  thick.  In  the 
XL200  stack,  201  bipolar  electrodes  were  aligned  back-to-back 
to  create  200  cells,  each  divided  by  a  membrane  0.2  mm  thick 
to  give  a  total  cell  thickness  of  about  7.5  mm.  The  Nation  mem¬ 
brane  material  was  assumed  to  have  the  properties  of  Nation  117 
(0.18  mm  thick,  and  with  a  conductivity  of  3Sm_1  [40]).  Each 
compartment  contained  a  HDPE  turbulence-promoting  mesh  to 
improve  mixing.  This  mesh  had  a  volumetric  porosity  of  about 
73%,  a  45-degree  cross  pattern,  a  strand  thickness  of  0.45  mm 
and  an  overall  thickness  of  about  0.90  mm.  The  XL200  stack 
was  bounded  by  two  mild  steel  endplates  coated  with  resin 
and  80  cm  wide,  155  cm  tall  and  4  cm  thick,  each  containing 
a  current  feeder.  Electrolytes  were  fed  to  the  cells  via  2.5-in. 
internal  diameter  chambers  formed  by  the  manifolds  of  the  elec¬ 
trode  frames,  and  through  spirals  in  the  manifolds  included  to 
reduce  shunt  currents.  The  electrolytes  were  then  returned  to 
large  storage  tanks  where  mixing  with  uncharged  electrolyte 
occurred.  A  heat  exchanger  was  included  on  the  less-corrosive 
sulphide  line,  which  maintained  the  electrolyte  temperature  at 
around  35  °C. 

4.2.  Flow  conditions  and  mass  transport 

The  variation  in  the  mass  transfer  coefficient  for  the  bromide 
compartment  of  an  XL5  module  (an  XL200  module  with  5  cells) 
has  been  found  to  be  [39]: 

Sh  =  0.081Re°'89Sc0'33  (22) 

for  20  <  Re  <120,  where  the  Sherwood  number  Sh  is  /<mde/D,  the 
Reynolds  number  Re  is  vde/v  and  the  Schmidt  number  Sc  is  v/D.  The 
hydraulic  diameter  de  was  1.9  mm.  The  relationship  for  the  slightly 
larger  sulphide  compartment  (valid  for  20  <  Re  <  150)  was  [41  ] : 

Sh  =  0.099Re°69Sc°33  (23) 

In  this  study,  the  average  electrolyte  velocity  v  was  fixed  at 
25mms-1  (a  typical  value  for  the  XL200  [19])  for  which  ReBr  and  Res 
were  62  and  46.  The  diffusion  coefficients  of  the  PSB  system  were 
estimated  to  be  12  x  10-6  cm2  s-1  (bromine),  5  x  10-6  cm2  s-1  (bro¬ 
mide),  6  x  10-6  cm2  s-1  (sulphide)  and  5  x  10-6  cm2  s-1  (sulphur) 
[20]. 

4.3.  Species  concentrations 

The  concentrations  of  the  oxidised  (S  and  Br2)  and  reduced 
(S2_  and  Br_)  were  calculated  on  each  side  of  the  membrane.  The 
effect  of  electro-osmosis  of  water  and  self-discharge  due  to  sulphide 
cross-over  were  included. 

Cations  crossed  the  cation-permeable  membrane  during  charge 
cycles  to  maintain  electro-neutrality,  dragging  with  them  associ¬ 
ated  water  molecules  that  resulted  in  significant  electrolyte  volume 
fluctuations.  Cations  crossed  from  the  bromide  to  the  sulphide 
compartment  on  charge  (decreasing  VBr  and  increasing  Vs),  and 
returned  on  discharge.  Sodium  ions  were  far  likelier  than  hydro¬ 
gen  ions  to  be  the  main  charge  carrier  in  the  PSB  system  due  to  its 
substantially  higher  concentrations  and  the  preference  of  Nation 
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membranes  for  sodium  ions  [40].  The  transference  coefficient  of 
water  tH2o  was  taken  as  10  molecules  of  water  per  cation  for  Nation 
117  carrying  pure  sodium  ions  [40]. 

Self-discharge  occurs  due  to  sulphide  ions  crossing  the  mem¬ 
brane  and  reacting  with  bromine  molecules: 

HS“  +  4Br2  +  4H20  ^  S042-  +  8Br“  +  9H+  (24) 

This  self-discharge  results  in  decay  of  the  equilibrium  potential,  and 
a  simple  estimate  of  the  decay  rate  was  by  differentiating  the  Nernst 
equation  for  the  PSB  system: 

dEeqm  RT  (  1  d[Br2]  1  d[S2“]  2  d[Br-]  1  d[S]\ 
“df- “2 F  ^IbT^]  dt  h[S5zr]  “ dt  [bF]- dt  [Si  ~dT J 

(25) 

The  rates  of  change  of  concentration  were  estimated  from  the 
fluxes  of  the  four  species  across  the  membrane  and  from  the  result- 


ing  reactions  that  occurred: 

d[Br2]  NA 
dt  =  VBr(  Js2~  hr2) 

(26) 

d[Br-]  NA 
dt  =Vb^+8JS2~  h'-] 

(27) 

d[S]  NA 
”dT  =  Vs  (“Js) 

(28) 

d  is2-]  j\M  /  1  \ 

(29) 

dt  =  VS  (  Js2-  4iBr0 

There  is  no  evidence  that  zero-valent  sulphur  crosses  the  mem¬ 
brane  in  significant  amounts.  The  diffusion  coefficient  of  molecular 
bromine  through  a  cation  exchange  membrane  with  a  water  sol- 
ventwasfoundtobe2.3  x  10-8  cm2  s-1  [42],  for  which  the  resulting 
bromine  flux  was  1.15  x  10-9  mol  cm-2  s-1  for  a  membrane  thick¬ 
ness  of  0.2  mm  and  a  bromine  concentration  of  1.0  M.  The  flux  of 
bromide  ions  across  a  Nafion  115  membrane  has  been  found  to  be 
3.6  x  10-9  mol  cm-2  s-1  [41].  However,  these  fluxes  were  consid¬ 
ered  to  be  negligible  compared  to  the  sulphide  flux.  Substituting 
Eqs.  (26)-(29)  into  Eq.  (25),  the  sulphide  flux  across  the  membrane 


was  found  from  the  voltage  decay  rate: 

_  _ 1000 _ 2 F 

S2‘  “  Wl((l/[S2-]rVs)  +  (4/[Br2 ]j-VBr)  +  (16/[Br-]rVBr))  & 


Experimental  studies  have  shown  that  the  long-term  decay  for 
an  XL200  stack  in  open-circuit  potential  can  be  taken  as  roughly 
constant  at  about  1.5  x  10-6  Vs-1  cell-1  [20].  Taking  the  open- 
circuit  voltage  to  be  equivalent  to  the  equilibrium  voltage,  the 
resulting  sulphide  flux  was  predicted  to  vary  significantly  during 
charge  with  an  expected  maximum  value  of  38  x  10-9  mol  cm-2  s-1 
(this  was  an  order  of  magnitude  larger  than  the  bromide  and 
bromine  fluxes  given  above,  confirming  that  the  sulphide  flux  was 
dominant). 

The  changes  in  species  concentration  along  the  electrode  were 
found  from  the  applied  current  and  sulphide  flux  as  follows: 


[Br2]h+i  [Br2lh+ 100QMchan  (±2  F  4/s2  ) 

(31) 

IBr-W,  -|Br-,h+  ,000^  (±J  +8V) 

(32) 

[S]h+1  -  [S]h  +  1000Mchan  (*2  F  ) 

(33) 

[S  ]h+l  -  [S  ]h+  -,0001/Achan  (±2  F  Js2~) 

(34) 

The  variation  of  the  tank  electrolyte  concentration  and  volume 
were  found  by  numerical  integration  of  the  differential  material 
balance  equations. 


4.4.  Numerical  optimisation  and  accuracy 

The  key  numerical  parameters  in  the  model  are  the  time  step 
At,  the  number  of  electrode  segments  H,  and  the  convergence  cri¬ 
teria  Ah  and  Xe.  Decreasing  the  size  of  time  step  At  had  a  4th-order 
effect  on  numerical  accuracy  (see  Fig.  3),  as  expected  since  Simp¬ 
son’s  rule  was  used  to  calculate  time-averaged  quantities.  At  low 
values  of  At  rounding  errors  prevented  further  improvement  in  the 
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Fig.  3.  Effect  of  numerical  parameters  on  estimated  relative  error  in  performance  measure:  ( x )  n,  (0)  Edis.f.  (□)  idis.entry.f  and  (A)  i'dis,exit,f-  Unless  specified  otherwise,  At=  1  min, 
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Fig.  4.  Effect  of  numerical  parameters  on  computational  effort:  (x)  Nit,  (+)  Tcomp-  Numerical  parameters  as  specified  in  Fig.  3. 


accuracy.  Increasing  the  number  of  electrode  segments  H  had  a  2nd- 
order  effect  on  accuracy,  while  the  current  and  voltage  convergence 
criteria  Ah  and  A.e  had  lst-order  effects. 

The  numerical  parameters  also  affected  computational  effort. 
Computational  cost  was  evaluated  by  recording  the  total  number 
of  iterations  required  for  convergence  and  the  computational  time 
required  to  run  the  model  in  Visual  Basic  on  a  Dell  Optiplex  GX620 
computer  with  a  Pentium  4,  3.2  GHz  processor  and  1024  Mb  of 
RAM.  As  can  be  seen  in  Fig.  4,  the  time  step  At  and  the  number 
of  segments  H  had  approximately  1  st-order  effects  on  the  number 
of  iterations  and  modelling  time  required.  Iteration  convergence 
criteria  had  little  effect  on  computational  effort  with  power  law 
exponents  for  Ah  and  Xe  of  only  0.03  and  0.05  respectively. 

This  information  is  summarised  in  Table  1.  Also  shown  are  the 
values  of  At,  H,  Xh  and  Xe  selected  on  the  basis  of  good  numerical 
accuracy  without  excessive  computational  time.  Numerical  errors 
in  integral  quantities  such  as  plant  efficiency  were  limited  to  10-4, 
with  the  error  for  segment  current  density  around  10-2. 

5.  Results  and  discussion 

The  model’s  results  were  compared  against  two  sets  of  exper¬ 
imental  data  for  an  XL200  stack  to  investigate  the  model’s  ability 
to  reproduce  observed  behaviour  under  different  operating  condi¬ 
tions. 


Table  1 

Orders  of  effect  of  numerical  parameters  on  model  accuracy  and  values  required  to 
limit  relative  numerical  error  to  10-4.  Numerical  parameters  as  specified  in  Fig.  3. 


At 

H 

A,h 

Ae 

Final  discharge  voltage 

0 

2nd 

1st 

1st 

Final  entry/exit 
discharge  current 

0 

1st 

1st 

1st 

density 

Energy  efficiency  n 

3rd-4th 

2nd 

1st 

1st 

Computational  time 

0.9 

1st 

0.03 

0.05 

Total  number  of 

0.9 

1st 

0.03 

0.05 

iterations 

Value  to  limit 
numerical  error  to 

1  x  10“4 

24  time  intervals 

30 

1  x  10“6V 

1  x  10~3  Am-2 

5.1.  Six-cycle  series  charged  at  600 A  mr2 

The  first  set  of  experimental  data  consisted  of  six  consecutive 
cycles  charged  for  65  min  and  discharged  for  60  min  at  a  cur¬ 
rent  density  of  600  Am-2.  During  each  changeover  the  stack  was 
held  under  open-circuit  conditions  for  10  min  to  allow  the  volt¬ 
ages  to  stabilise.  A  summary  of  the  parameters  used  to  model  this 
series  of  XL200  PSB  cycles  is  given  in  Table  2.  A  detailed  discus¬ 
sion  of  the  selection  of  these  parameters  can  be  found  elsewhere 
[20]. 

Fig.  5  shows  the  predicted  voltage  profile  obtained  from  the 
model,  along  with  the  experimental  voltage  trace  recorded  for  a 
working  XL200  PSB  system.  The  model  results  agreed  well  with  the 
experimental  data  in  a  number  of  ways.  In  particular,  the  average 
value  of  the  energy  efficiency  was  calculated  as  59.14  and  59.16% 
for  the  model  and  experiment  respectively  [20].  Both  the  actual 
and  predicted  open-circuit  voltages  decayed  over  time  at  similar 
rates  due  to  self-discharge,  and  by  the  end  of  the  series  the  sys¬ 
tem  was  very  near  to  mass-transport  limits  as  can  be  seen  by  the 


Fig.  5.  Predicted  and  recorded  cell  voltage  for  a  six-cycle  series  of  an  XL200  PSB 
system  with  the  input  data  given  in  Table  2:  (a)  predicted  voltage  and  (b)  recorded 
voltage. 
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Battery  design  used  for  modelling  a  series  of  6  consecutive  PSB  cycles  charged  at  400  Am-2. 
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System  properties 

Symbol  (unit) 

Value 

Electrolyte  velocity 

v(mms-1) 

25 

Number  of  cells 

N  (cells) 

200 

Transference  coefficient 

t 

10  molecules  H20  per  Na+  ion  transferred 

OCV  voltage  decay  rate 

(Vs-1  cell- 

5  x  10-7 

Half-cycle  properties 

Charge 

Discharge 

Applied  current  density 

ich/idis  (Am-2) 

600 

600 

Cycle  length 

tch/tdis  (min) 

65 

60 

Half-cell  properties 

System  A 

System  B 

Oxidised/reduced  on  charge 

- 

Oxidised 

Reduced 

Standard  redox  potential 

E°  (V) 

+1.09 

-0.48 

Kinetic  rate  constant 

ks  (xlO-8  ms-1) 

40 

3 

Initial  electrolyte  volume 

V0  (m3) 

3.5 

3.5 

Electrolyte  viscosity 

v  (xlO-6  m2  s_1) 

0.76 

1.04 

Initial  state  of  charge 

0.3 

0.19 

Mass  transport  relationship 

- 

Sh  =  0.081  Re089  Sc0,333 

Sh  =  0.099Re°-69Sc°-333 

Species  properties 

Bromine 

Bromide 

Sulphur 

Sulphide 

Stoichiometric  coefficient 

z 

1 

2 

1 

1 

Diffusion  coefficient 

D  (xlO-10  m2  s-1) 

20 

12 

5 

6 

Initial  electrolyte  concentrations 

[P]  (M) 

0.73 

3.42 

2.86 

1.53 

Electrolyte  properties 

System  A 

System  B 

Membrane 

Conductivity 

k  (Snrr1) 

23 

29 

1 

Thickness 

l  (mm) 

0.95 

0.95 

0.2 

sharp  drop-off  in  cell  voltage  towards  the  end  of  discharge.  There 
were  also  a  number  of  differences  in  the  data.  The  predicted  voltage 
changed  more  slowly  during  half-cycles  (especially  on  charge).  This 
could  be  because  side  reactions  such  as  gas  evolution  were  ignored 
and  the  model  could  not  exceed  mass-transport  limits.  The  model 
over-predicted  cell  voltage  on  charge,  and  this  is  probably  due  to 
the  use  of  the  same  rate  constant  for  both  forward  and  reverse 
reactions.  The  model  open-circuit  voltage  stabilised  immediately 
as  capacitance  effects  and  the  time  taken  to  reach  equilibrium  were 
neglected. 

The  predicted  variation  in  species  concentration  in  the  elec¬ 
trolyte  tanks  is  shown  in  Fig.  6.  The  rate  of  change  of  concentration 
varied  during  half-cycles  due  to  the  varying  cross-membrane  sul¬ 
phide  flux  and  the  change  in  electrolyte  volume  (due  to  the 
electro-osmotic  water  flux).  The  concentration  of  Br2  decreased  and 
the  Br-  concentration  increased  during  the  series  of  cycles  due  to 


Time  (hr) 

Fig.  6.  Predicted  variation  in  electrolyte  tank  concentrations  for  a  six-cycle  series 
of  an  XL200  PSB  system  with  the  input  data  given  in  Table  2:  (a)  [Br2],  (b)  [Br~  ],  (c) 
[S2_]  and  (d)  [S]. 


a  combination  of  self-discharge  and  the  net  transfer  of  water  to 
the  sulphide  electrolyte  as  the  duration  of  the  charge  cycle  tc h  was 
longer  than  the  discharge  cycle  tdis.  Over  the  series  of  cycles,  the 
concentration  of  S2-  fell  slightly  due  to  the  sulphide  flux  through 
the  membrane,  and  the  zero-valent  sulphur  concentration  also  fell 
slightly  due  to  the  longer  charge  period. 

The  resulting  variation  in  the  state  of  charge  is  shown  in  Fig.  7. 
Both  cr s  and  <rB r  increased  on  charge  and  decreased  on  discharge, 
but  neither  electrolyte  used  its  full  state  of  charge  range.  The  value 
c rBr  was  limited  in  order  to  minimise  the  risk  of  corrosion  and  Br2 
release.  The  initial  value  of  (jBr  was  selected  to  ensure  that  it  did 
not  decrease  to  zero  due  to  self-discharge.  Greater  utilisation  of  the 
sulphide  electrolyte  could  be  achieved  by  reducing  the  electrolyte 
volume;  the  resulting  higher  species  concentrations  would  improve 
mass  transport,  and  smaller  electrolyte  and  tank  costs  would  be 
required.  The  fluctuation  of  crs  during  cycles  was  less  than  that  of 
cr Br  due  to  the  greater  amount  of  total  sulphur  in  solution.  Even 


Fig.  7.  Predicted  variation  in  state  of  charge  for  a  six-cycle  series  of  an  XL200  PSB 
system  with  the  input  data  given  in  Table  2:  (a)  crBr  and  (b)  <ts. 
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though  the  duration  of  the  charge  cycle  tch  exceeded  the  discharge 
cycle  tdis,  the  bromide  state  of  charge  o'er  decreased  with  each  cycle 
due  to  the  sulphide  flux.  Conversely  the  sulphide  state  of  charge  crs 
rose  slightly  over  time  due  to  a  decrease  in  the  concentration  of 
zero-valent  sulphur. 

The  predicted  variation  in  the  current  density  is  shown  in  Fig.  8 
at  a  number  of  locations  along  the  electrode.  The  current  density 
was  lower  at  the  electrode  outlet  than  at  the  inlet  due  to  the  lower 
concentrations  of  reactant.  The  range  of  current  density  variation 
along  the  electrode  increased  during  half-cycles  as  the  conversion 
per  pass  increased,  increasing  the  mass  transport  overpotential  and 
hence  reducing  the  local  current  density  at  the  outlet.  The  range  of 
current  densities  was  also  larger  on  discharge  than  on  charge  as 
the  system  was  closer  to  mass  transport  limiting  conditions.  This 
range  of  current  densities  on  discharge  increased  each  cycle  as  the 
bromide  side  state  of  charge  decreased  leading  to  increasing  mass 
transport  overpotentials  at  the  cell  outlet. 

The  predicted  variation  in  the  dimensionless  current  density,  id 
(the  ratio  of  current  density  to  limiting  current  density),  is  shown 
in  Fig.  9.  As  expected,  the  current  density  at  the  exit  is  closer  to 
mass  transport  limiting  conditions  than  at  the  inlet.  During  each 
half-cycle,  id  increased  as  the  limiting  current  decreased  with  reac¬ 
tant  concentrations.  For  most  of  the  six-cycles,  the  current  density 
was  always  less  than  50%  of  the  limiting  current.  However,  dur¬ 
ing  the  last  few  cycles,  limiting  conditions  are  approached  at  the 
end  of  discharge  on  the  bromide  side.  This  led  to  high  overpoten¬ 
tials  and  reduced  efficiency,  and  consequently  the  experiment  was 
terminated. 

5.2.  Seventeen-cycle  series  charged  at  400 A  m~2 

A  series  of  17  consecutive  charge-discharge  cycles  was  also  stud¬ 
ied.  This  series  began  with  3.8  and  3.2  m3  of  fresh  bromide  and 
sulphide  electrolytes  respectively,  and  the  initial  species  concentra¬ 
tions  were  4.5  M  Br-,  3.8  M  S  and  1  M  S2-.  For  modelling  purposes, 
the  initial  Br2  concentration  was  assumed  to  be  10-4  M  as  zero  con¬ 
centrations  are  incompatible  with  the  Nernst  equation.  The  charge 
and  discharge  current  densities  were  400  and  600 Am-2  respec¬ 
tively.  For  the  first  cycle  the  charging  time  was  175  min,  and  for 
subsequent  cycles  it  was  110  min.  The  discharge  time  was  65  min  for 
each  cycle  except  for  the  twelfth  cycle,  for  which  58.5  min  was  used 
for  balancing  purposes.  Between  each  changeover  the  stack  was 
held  under  open-circuit  conditions  for  10  min  to  allow  the  voltage 
to  stabilise. 


Time  (hr) 

Fig.  9.  Predicted  variation  in  dimensionless  current  id  =  i/i'l  for  a  six-cycle  series  of 
an  XL200  PSB  system  with  the  input  data  given  in  Table  2:  (a)  id,Br,entry.  (b)  id.Br.exit- 

(c)  id,S, entry  ^nd  (d)  l'd,S,exit- 

For  the  model,  lower  rate  constants  /<s  Br  and  ks  S  of  10-7  and 
10-8ms-1  respectively  were  selected  in  order  to  avoid  disconti¬ 
nuities  in  overpotential  occurring  at  the  transition  point  between 
linear  and  Tafel  approximations  to  the  Butler-Volmer  equation  [20]. 
It  was  not  possible  to  achieve  a  current  density  on  discharge  of 
600 Am-2  without  exceeding  limiting  current  conditions.  It  was 
therefore  assumed  that  a  constant  current  efficiency  of  93.3%  was 
achieved  (i.e.  a  discharge  current  density  of  560  Am-2  was  applied 
in  the  model).  A  lower  OCV  decay  rate  of  4  x  10-7  V  s-1  cell-1  was 
used  in  order  to  prevent  species  concentrations  falling  to  zero  over 
the  17-cycle  series.  A  membrane  conductivity  of  3  S  m-1  was  used; 
other  model  settings  were  as  given  in  Table  2. 

The  resulting  model  cell  voltage  is  shown  in  Fig.  10  along  with 
the  experimental  cell  voltage  data  for  comparison.  The  model  over¬ 
estimated  the  voltage  during  charging,  probably  due  to  the  low 
rate  constant  used.  Conversely  on  discharge  the  cell  voltage  was 
underestimated  due  to  the  lower  discharge  current  density.  The 
predicted  rate  of  change  of  voltage  during  half-cycles  was  also 
lower  than  the  observed  rate.  Nevertheless  the  agreement  was  good 
in  several  aspects.  In  particular  the  sharp  drop-off  in  cell  voltage 
on  discharge  was  observed  in  both  the  model  and  experimental 
data  due  to  mass  transport  overpotentials.  The  experimental  data 
indicated  that  higher  mass  transport  overpotentials  occurred  than 
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Fig.  8.  Predicted  variation  in  current  density  distribution  for  a  six-cycle  series  of  an 
XL200  PSB  system  with  the  input  data  given  in  Table  2:  (a)  entry  segment,  (b)  5th, 
(c)  10th,  (d)  15th,  (e)  20th,  (f)  25th  and  (g)  exit  segment. 


Fig.  10.  Predicted  and  recorded  variation  in  cell  voltage  for  a  17-cycle  series  of  an 
XL200  PSB  with  the  input  data  specified  in  Section  4.2:  (a)  predicted  voltage  and  (b) 
recorded  voltage. 
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Fig.  11.  Predicted  variation  in  state  of  charge  for  a  17-cycle  series  of  an  XL200  PSB 
with  the  input  data  specified  in  Section  4.2:  (a)  aB r  and  (b)  crs. 

those  predicted  by  the  model.  Simulations  of  the  cell  voltages  close 
to  mass  transport  limiting  conditions  were  found  to  be  very  sensi¬ 
tive  to  the  applied  current.  The  shorter  discharge  period  on  cycle  12 
was  used  to  increase  the  average  state  of  charge  and  hence  reduce 
mass-transport  limits  overpotentials  during  discharge.  Although 
the  effect  of  this  shorter  cycle  in  later  cycles  is  evident  in  the  model 
results,  it  is  less  clear  in  the  experimental  data. 

The  predicted  variation  in  state  of  charge  is  shown  in  Fig.  11.  On 
the  bromide  side,  txBr  cycled  roughly  between  0.1  and  0.3  through¬ 
out  the  series  as  the  loss  of  bromine  due  to  the  sulphide  flux  was 
counterbalanced  by  the  net  excess  of  bromine  resulting  from  the 
greater  amount  of  charge  stored  on  charge  than  was  recovered  on 
discharge.  The  low  state  of  charge  on  the  bromide  side  is  associated 
with  the  high  mass  transport  overpotentials  at  the  end  of  discharge 
as  the  Br2  concentration  falls.  On  the  sulphide  side,  crs  steadily 
increased  during  the  series  due  to  the  decay  in  [S].  Towards  the 
end  of  the  12  cycles,  mass  transport  overpotentials  are  evident  at 
the  end  of  charge,  due  to  decreasing  concentrations  of  zero-valent 
sulphur. 

The  predicted  variation  in  id  is  shown  in  Fig.  12.  On  the  bromide 
side,  id  Br  was  very  low  during  charge  throughout  the  series  due 
to  the  high  concentrations  of  bromide.  On  discharge  id >Br  was  very 
high,  frequently  exceeding  0.9  at  the  electrode  exit.  This  confirms 
that  the  rapidly  decreasing  potentials  at  the  end  of  discharge  are  due 
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Fig.  12.  Predicted  variation  in  dimensionless  current  id  for  a  17-cycle  series  of  an 
XL200  PSB  with  the  input  data  specified  in  Section  4.2:  (a)  i'd,Br, entry.  (b)  i d,Br,exit.  (0 

id, S, entry  ^fld  (d)  id.S.exit- 


to  increasing  mass  transport  overpotentials  at  the  bromide  elec¬ 
trode.  After  cycle  12  id  Br  rarely  exceeded  0.6  as  the  shorter  discharge 
period  moved  the  system  away  from  mass-transport  limits. 

On  the  sulphide  side,  ids  was  initially  higher  on  discharge 
than  on  charge  due  to  relatively  low  sulphide  concentrations.  Over 
time,  id  s  became  higher  on  charge  due  to  falling  concentrations 
of  zero-valent  sulphur.  Again,  this  confirms  that  mass  transport 
overpotentials  occur  during  charging  in  later  cycles  at  the  sulphide 
electrode. 

6.  Conclusions 

This  study  has  shown  that  modelling  of  redox  flow  battery  sys¬ 
tems  for  energy  storage  can  reveal  important  characteristics  and 
loss  mechanisms  over  a  series  of  cycles  under  different  operating 
conditions.  For  example,  for  the  conditions  studied,  mass  trans¬ 
port  overpotentials  at  the  bromide  electrode  were  found  to  limit 
the  performance  during  discharge.  Enhancement  of  mass  trans¬ 
port  coefficients  in  the  bromide  electrolyte  may  be  recommended. 
The  model  showed  that  significant  drift  in  conditions  could  occur 
due  to  self-disharge  and  electro-osmotic  effects.  Careful  electrolyte 
management  will  be  required  for  reliable  operation  of  the  PSB  RFB 
system.  The  model  can  also  be  used  to  find  suitable  electrochem¬ 
ical  rate  constants  by  fitting  with  experimental  data.  For  the  PSB 
system  at  the  composite  activated  carbon  electrodes  developed  by 
Regenesys,  bromide  and  sulphide  rate  constants  of  4  x  10-7  and 
3  x  10-8ms-1  were  obtained.  In  a  following  paper  [18]  we  show 
how  the  model  can  be  used  for  design  and  optimisation  of  a  large 
scale  RFB  system. 
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